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homogeneous and isotropic flat FLRW model. By proper choice of the particle creation 
rate as a function of the Hubble parameter, the model shows the evolution of the universe 
starting from the inflationary scenario to the present accelerating phase, considering the 
cosmic matter as normal perfect fluid with barotropic equation of state. 


Keywords: Cosmic evolution. Particle creation, Bulk viscosity, Isentropic process. 

Cosmography 

PACS Numbers: 04.70.Dy, 98.80.-k, 04.90.-he 


^ schakraborty@math.jdvu.ac.in 
^ span@research.jdvu.ac.in 
subhajitl729@gmail.com 



2 


I. INTRODUCTION 


In cosmology, usually homogeneous and isotropic flat Friedmann-Lemaitre-Robertson-Walker 
(FLRW) model is considered due to its agreement with inflation and the cosmic microwave back¬ 
ground (CMB) observations. Here, the only dissipative phenomenon is in the form of bulk viscous 
pressure which may occur either due to coupling of different components of the cosmic substratum 

nQ 

iHSj, or, due to non-conservation of (quantum) particle number [GHaj. In the present work, second op¬ 
tion is only considered and for simplicity of calculations, isentropic (i.e., adiabatic) particle production 
IQ of perfect fluid particles is considered, and, consequently, there is a simple linear relationship 
between particle production rate, and, the viscous pressure. However, it is to be noted that still there 
is entropy production due to enlargement of the phase space (due to increase in the number of fluid 
particles, and, also, due to the expansion of the universe in the present model) of the system. 

Further, in the context of non-equilibrium thermodynamical prescription, second order deviations 
from equilibrium (following Israel and Stewart ll|, [l^ ) are considered in order to eliminate the draw¬ 
backs of the 1st order Eckart theory related to causality and stability. As a result, entropy flow will be 
related to the cosmological particle production of isentropic nature, and, the bulk viscous pressure will 
become a dynamical degree of freedom having causal evolution equation. Furthermore, it is possible 
to eliminate bulk viscous inflation without particle production (as used by many authors), and, thus, 
one may have a model universe starting from a de Sitter phase which gradually evolves to standard 
FLRW model. The present work will be an attempt to incorporate the present accelerating phase 
within this causal second order theory. 

Usually, to explain the recent accelerated expansion of the Universe as predicted by Supernovae 
la and complementary observations Q, 14], there are two common approaches — one within the 
framework of Einstein gravity by introducing an unknown type of matter component with a large 
negative pressure (dark energy), and, secondly, by modifying Einstein gravity theory, and interpreting 
the extra geometric terms as hypothetical matter component to explain the present accelerating phase. 
However, both the approaches are concentrated only to explain the recent observational predictions 
— there is no concern about the past or future evolution of the universe. In this context, the present 
work is an attempt not only to explain the recent observations, but also the past evolution of the 
universe without following any one of the above mentioned conventional approaches. 


The paper is organized as follows: Section |TI] deals with a brief review of the non-equilibrium ther¬ 
modynamics in cosmology, section IHII describes the explicit solutions of the cosmological parameters 
corresponding to a unique particle creation rate. In section HYi we have analyzed the cosmographic 
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parameters associated with our model. Section 0 we have described the particle productions in the 
language of field theory. In section m we present an associated Hawking like radiation. Finally, we 
finish our discussions in section IVIII 


II. PARTICLE CREATION MECHANISM AND NON-EQUILIBRIUM 
THERMODYNAMICS: BASIC EQUATIONS 

In a closed thermodynamical system, suppose there are N particles having internal energy E. Then 
the first law of thermodynamics which is essentially the conservation of internal energy reads 


dE = dQ-pdV, (1) 

where dQ is the amount of heat received by the system in time dt, and as usual p is the thermody¬ 
namic pressure, V is any co-moving volume. The above energy conservation relation can be rewritten 
as Gibb’s equation 


Tds = dq = d 


+ pd 


(2) 


where ‘s’ is the entropy per particle (specific entropy), T is the temperature of the fluid, n = N/V 
is the particle number density, and dq = dQ/N is the heat per unit particle. It should be noted that 
the above Gibb’s equation also holds for open thermodynamical system, i.e., when the particle number 
is not conserved [o, 2, 

As mentioned earlier, we consider spatially flat FLRW model of the universe as an open ther¬ 
modynamical system, and, non-equilibrium nature comes into picture due to the particle creation 
mechanism. Thus, the Einstein’s field equations take the form 


3H^ = up, and, 2H = —k{p + p+ 11), (3) 

where k = 8ttG is the Einstein’s gravitational constant. Eor a relativistic fluid with dissipation in 
the form of bulk viscosity has the energy-momentum tensor 


= {p + p + + {p + 


( 4 ) 


and the conservation equation = 0 (Bianchi’s identity) reads 
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p + 6{p + p + Il) = 0, (5) 

where the bulk viscous pressure 11 is related to the entropy production. 

Further, in an open thermodynamical system, the non-conservation of fluid particles is reflected by 
the equation 


= h + 9n = nV , 


( 6 ) 


where F denotes the rate of change of the particle number (A^ = no?) in a co-moving volume a?. 
= nu^ is the particle flow vector, n is the particle number density, is the unit time-like vector 
(4—velocity), 6 = u^-f^ is the expansion of the congruence of time-like geodesics and notationally 
h = naU^". So, creation and annihilation of particles is characterized by the sign of F (creation: 


F > 0 and annihi 


pressure 


ation: F < 0). Further, a non-zero F is dynamically equivalent to an effective bulk 


16 h 21 || of the fluid, and, hence, non-equilibrium thermodynamics comes into consideration. 


However, Lima et al. 


22l | showed that such scalar processes (bulk viscosity and matter creation) are not 


equivalent from a thermodynamic viewpoint — only the dynamic behavior can simply be demonstrated 
in the case of “adiabatic” particle creation as follows: 

Using the above conservation equations ([5]) and ([6]) into the Gibb’s relation (i.e., Eq. ([2])), one 
obtains the entropy variation per particle as 


( 7 ) 

For simplicity, if the thermal process is assumed to be isentropic (i.e., adiabatic) then the entropy 
per particle remains constant (in contrast to dissipative process), i.e., s = 0 and we have 

n = -^(p + p)- (8) 

The above relation shows that a non-vanishing F will produce an effective bulk pressure on the 
thermodynamic fluid and non-equilibrium thermodynamics comes into the picture. In other words, 
a dissipative fluid is equivalent to a perfect fluid having a non-conserved particle number. It should 
be noted that there is entropy production only due to the enlargement of the phase space of the 
system. Further, one may note that, this effective bulk pressure does not correspond to conventional 


s = — 


nT 


^ + -q{p + p) 
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non-equilibrium phase, rather a state having equilibrium properties as well (but not the equilibrium 
era with F = 0). 


In the framework of second order non-equilibrium thermodynamics due to Israel and Stewart 
the entropy flow vector (S"^) has the expression 


Q, 


rll^ 

= sIV“ - —n“, (9) 

2cr ’ ^ ^ 

where C is the coefficient of bulk viscosity and r represents the time of relaxation. Now, using the 
conservation equations ([5]) and ([U]), one obtains the entropy production density from the above Eq. 
(HI) as 


TS^-a = -nfiV - n 


tFI 

H-h -iir 

C 2 



;a 


( 10 ) 


where /j = 
thermodynamics (i.e., 


— Ts is the chemical potential. Now, for the validity of the second law of 
5'“;a > 0) one can choose the ansatz for bulk viscous pressure as 


n = -c 


r • 1 

' + -11 + -nr 

c 2 


CT 


fj,nT 


Cn 


( 11 ) 


As a result, 11 becomes a dynamical variable with an inhomogeneous evolution equation 


n 2 + rnn + + (ue = ( 12 ) 

Here, the chemical potential /i may act as an effective symmetry-breaking parameter in relativistic 
field theories. Also, the evolution equation becomes linear hrst order in nature in absence of chemical 
potential. 

In this context, it is relevant to mention that the basic physical difference between the noncausal 
and the causal theory is the introduction of the time of relaxation (in the later one). As a result, in 
causal theory H decays to zero after F has been switched off (assuming non-vanishing F produces the 
effective viscous pressure). Moreover, if the above second order theory is isentropic in nature then the 
entropy production density simplihes to 


= -n 


[n 

f T \ 

r ■ 

nsT 

[2 


H—n - 

;a CT 

n 


(13) 












6 


Further, due to the isentropic condition ([5]), the evolution of the relevant thermodynamical variables 


are 


23| 


p =-[e-T){p + p), p =-cl{e-T){p + p), - = -{e-T), = (u) 

n lop 


where ^ is the square of the adiabatic sound velocity 23|. 


III. PHENOMENOLOGICAL CHOICE OF PARTICLE CREATION RATE AND COSMIC 

EVOLUTION 

To describe the cosmic evolution (for a given particle creation rate as a function of the Hubble 
parameter), one can eliminate the effective bulk pressure H from the Einstein field equations ([3]) using 
the isentropic condition ([8]) to obtain 


T = i + A 4,' 

3H 37 I 


(15) 


On the other hand, this equation can also be considered as the determining equation for the particle 
creation rate F from the given cosmic evolution. Also, the deceleration parameter q takes the form 




(16) 


So, in the present context, the cosmic history is characterized by the fundamental physical quanti¬ 
ties, namely, the expansion rate H and the energy density /?, and, as a result, the gravitational creation 
rate F can be defined in a natural way. 

In earlier studies 2214271] . the particle creation rate F has been chosen for different phases of the 


evolution from thermodynamical viewpoint. As F should be greater than H in the very early universe 
so that the created radiation may be considered as a thermalized heat bath. Hence, at the very early 


Universe HT is c 
solution [22, 23 


losen to be proportional to 


m, 


(i.e., 


F oc p) and the corresponding cosmological 


27( 1 shows a smooth transition from the inflationary scenario to the radiation era. 


Also for this adiabatic production of relativistic particles, the energy density scales as pr oc T^, i.e.. 


black body radiation 


23 , 12711 . 


Furthermore, it has been recently shown 


2S 


29(1 that, T (X H and T x 1/H describe respectively 


the intermediate matter dominated era (starting from radiation) and the transition from matter dom¬ 
inated era to late time acceleration. The scale factor, the Hubble parameter and the thermodynamical 
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FIG. 1. This figure describes the scale factor over the cosmic time. 



FIG. 2. The figure shows the behavior of the 
Hubble parameter over the cosmic time. 



FIG. 3. This displays the entire cosmic history 
from inflation to present accelerating phase of the 
universe. 


parameters are shown to be continuous across the transition points (i.e., the epochs from the infla¬ 
tionary era to the radiation era and from the matter dominated era to the late time acceleration). 


Also, r = Fq has been shown m to correspond the emergent scenario. Subsequently, in another work 
a linear combinations of all these choices, i.e., F = Fq + IH^ + mH + n/H has been examined 


to describe the whole evolution of the Flniverse. Although no exact analytic solution is possible for 
this choice of F (from Eq. (Ilbjl l. still the graphical representation of the deceleration parameter has 
shown the whole evolution of the universe starting from an early inflationary epoch to the present 
accelerating phase, and, the model predicts a possible transition from present accelerating stage to 
decelerating phase again in the future. 
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With this background in mind, the present work is a partial modification of the above general choice 
with an aim to have an exact analytic solution so that more sophisticated cosmic study, namely, the 
cosmographic analysis can be done. By choosing the coefficients appropriately (which will be clear 
subsequently) the form of T is taken as 


r 


n a 

—/r + 3H + — 


(17) 


with /i and a as real constants. Now substituting this T into the evolution equation (I15p . the 
explicit solution reads as 


a = aoe 




exp 




H = 


a 




H— 7 ^ exp 






(18) 

(19) 


where oq and to ^.re the constants of integration. The graphical representation of the cosmic 
evolution, namely, the scale factor a, the Hubble parameter H, and the deceleration parameter q are 
presented in FIGs. (1)~(3) respectively for various choices of 7 , the equation of state parameter for the 
cosmic fluid. The diagramatic representation of q shows two transitions of q (from accelerating phase 
to deceleration, and, then again acceleration) which indicates that the present model of the Universe 
describes the evolution from the inflationary scenario to the present late time acceleration through 
the decelerated matter dominated era. Thus, we have a complete cosmic history after the big bang 
till today. It is interesting to note that as the cosmic time becomes very large, we have the ACDM 
model: 


a ~ aoe'^°*, 77 ~ Ffg = 1 Q — —1; j P — 377o = A = —p, 


( 20 ) 


and it agrees with the recent Planck data set 


32l |. Now, corresponding to the cosmological solution 


(i.e., Eqns. m and (jl 8 p i of the present model, the relevant thermodynamical parameters evolve as 
(see Eq. IfTTp i 


r = rn 


n = no 


a i 
—^ H — T^exp ■ 
p^ p^ 


1 

—^ 7 —77 exp 

p^ p^ 




(t - to) 


(^) 


-^(t-to) 


S 


Soe 



exp 





a i 
^ + -^exp 

p^ p^ 




( 21 ) 
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FIG. 4. This is the behavior of the temperature 
over the cosmic time for three different choices 
of 7 . 

where So is some constant, and, Tq, tiq are integration constants. We have plotted these thermo¬ 
dynamical parameters in FIGs. (4)-(6) for various choices of 7 . It should be mentioned that if 7 = 2, 
i.e., p = p (relativistic fluid), then T oc pi, which represents the usual black body radiation. Note 
that, in the ACDM limit the above thermodynamical parameters become 



which shows that although the temperature and the number density become constant, but the 
entropy in a co-moving volume evolves as the scale factor. Lastly, we note that, in the above limit, 
the particle creation rate becomes constant: F —)• Fq = 2>Hq. 


IV. COSMOGRAPHIC ANALYSIS AND OBSERVATIONAL DATA 


In this section, we make a comparative study of the present model of the universe with the presently 


available observed data set. At first, we discuss a geometric view of the DE models. Sahni et ah 


M 


first proposed this idea with two dimensionless and model independent geometric parameters {r, s} 
defined as 


1 d^a 
aH^ dt^ ’ 



r = 


and, s 


(23) 
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FIG. 5. The figure shows the variation of the 
number density with the evolution of the 
universe. 



FIG. 6. This is the entropy variation of the 
universe with its evolution. 


These two parameters in Eq. (1231) are used to filter the observationally supported DE models 
from other phenomenological DE models existing in the literature. Subsequently, these geometric 
investigation was further extended by considering the Taylor series expansion of the scale factor about 
the present time as 


^ = 1 + - tp) + - hf + 0[{t - t,)\ (24) 

where the model independent parameters j, s, I, m are known as cosmographic parameters [^, 
and, are defined as 


J = 


1 cfa 

am dt^ ’ 


s = 


1 

aH^ dt^ ’ 


I = 


1 d^a 

am dt^ 


and, m = 


1 d^a 

am dt^ ' 


(25) 


Here, the suffix ‘p’ stands for the value of the corresponding variable at the present epoch (tp). It 
should be noted that, the cos mog raphic parameters are individually named as jerk (j) (this ‘j’ is same 
as ‘r’ dehned by Sahni et a l. |33l] b snap (s) (this ‘s’ is different from one defined by Sahni et al. 33|), 
lerk, and m parameter j^, Ssl]. 

The above cosmographic parameters (CP) can be expressed in terms of the deceleration parameter 
(g), and, its higher derivatives in the following way: 




(26) 
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(27) 


(28) 


(29) 



FIG. 7. The figure shows a comparative study of 
the deceleration parameter of our model with the 
4 different latest observational data sets. 

Further, the deceleration parameter can be expressed in terms of the redshift parameter 2 : (= 
1/0 — 1) as 


q{z) =qp + i-qp - 2q‘^ + jp)z + -{2qp + + 8q^ - 7qpjp - Ajp - Sp)z^ + 0(2^). (30) 


For the observed data sets, we choose the following 


(i) 192 Sne la and 69 GRBs with CPL parametrization (data 1) 361 

(ii) 192 Sne la and 69 GRBs with linear parametrization (data 2 ) |36 | 
(hi) Supernovae Union 2+ BAO+ OHD+ GRBs data (data 3) 37] 
(iv) Supernovae Union 2+ BAO+ GRBs data (data 4) 


In FIG. 7, we have made a comparative study of the deceleration parameter for the present model 
with those for the above four data sets. From the graph we see that the behavior of the deceleration 
parameter in our model almost matches with the 4 latest observed data sets. 
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t 

- y = 0.9-Y = 1 “ ■ “ y = 1.1 


FIG. 8 . The figures show the variation of 
the “jerk” (j) throughout the evolution of 
the universe. 



FIG. 9. This is just to show that during the 
early phase of the universe, ‘j’ started from 
— ve values. 


Also, we have shown the graphical representations of the above four CP parameters for different 
choices of 7 . 

In fact, FIG. 8 shows the complete evolution of the jerk parameter ‘j’, while FIG. 9 shows explicitly 
the variation of j in the early phase of the universe. There is a transition of j from — ve values to 
+ ve values. The variation of the snap parameter ‘s’ over the entire evolution of the universe has 
been shown in FIG. 10, while FIGs. 11 and 12 represent the detailed variation of ‘s’ in the early 
and late phases respectively. There are three transitions of ‘s’ during the whole evolution. FIG. 13 
shows the variation of the lerk parameter ‘F over the entire cosmic time, and its variation at early 
phase is presented in FIG. 14. The figures show two transitions of the lerk parameter. Finally, the 
‘m’ parameter is graphically represented in FIGs. 15-17, which also has the two transitions over the 
entire evolution. 


V. FIELD THEORETIC DESCRIPTION OF COSMIC HISTORY 

In this section, we shall describe the cosmic evolution from the field theoretic point of view by 
describing the whole dynamical process as the evolution of a scalar field (p having self interacting 
potential V{4>), or, equivalently, the evolution of the present effective imperfect fluid can be described 
by a minimally coupled scalar field. Thus, the energy density and the thermodynamic pressure of the 
cosmic fluid are given by 
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FIG. 11. The figures are just to show that 

FIG. 10. The figures show the variation of the , . , , , r i 

during the early phase of the universe, 

“snap” (s) parameter throughout the entire 

s started from — ve values. 

evolution of the universe. 



FIG. 12. This is just to show that when t —> oo, 
‘s’ becomes + ve. 


+ y{4>), and, = p + U = -V{4>). 

Thus, for the present isentropic thermodynamical system, we have 

= —{p?H — a^), 

K 


(31) 


(32) 

(33) 
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FIG. 13. The figures show the variation of 
the “lerk” { 1 ) parameter throughout the entire 
evolution of the universe. 


FIG. 14. The figures are just to show that 
during the early phase of the universe, ‘F 
was + ve. 


where 11 is eliminated by the isentropic condition equation ([ 8 ]), particle creation rate T is obtained 
from Eq. (HID, and the first Friedmann equation in ([3]) has been used. Thus, integrating the above 
equation ([32]) . the explicit form of cj) is 


(t) = (t>o -(34) 

and the potential can be expressed explicitly in (j) as 

V{cl>) = Vo + Vi{cl> - cl>of + V2{(P - 0o)^ (35) 

where cpQ is the constant of integration, and, Vq, Vi, V 2 are the constants depending on n, a, and 
7 . Note that, the scalar field has always a value less than (pQ. Further, the particle creation rate F 
can be expressed as a function of (j) as 


r - Fo + ri(</> - </>o)' + 

with Fq, Fi, and F 2 as constants. 

It is worthwhile to mention that, in the asymptotic limit (i.e., ACDM era), (j) becomes a constant 
(= (po), and the potential behaves as the cosmological constant A. 
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FIG. 15. The figures show the variation of the during the early phase of the universe, ‘to’ 
‘to’ parameter throughout the entire evolution started from + ve 


of the universe. 



- Y = 0.9 Y = 1 — “ Y " 1-1 


FIG. 17. This is just to show that when t —> oo, 
‘to’ becomes + ve. 


Furthermore, in view of the slow roll approximation during the inflationary phase, the density 


fluctuations are of the form Sh r-^ ~ 10 


-5 


1^, 27]. So, 


for the present model, we see 


^ kH . 

~ 7 „ \ llnitial epoch-. H=Hj j i.e., 

(l - ^) 


(37) 


(38) 
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As K = 8 ttG = Sirlpi (in units = c = 1), so, the ratio of the Hubble parameter at the two 


accelerating phases (initial and present) is proportional to . At the Planck size of the universe, 
i.e., Ipi = 10“^^ m, the Hubble parameter is Hj ~ 10^^ sec“^ SS], so, Hq ~ 10 sec Thus, the 


cosmological parameter A is ~ 10 the present observed value. 


VI. HAWKING LIKE RADIATION IN THE CONTEXT OF PARTICLE CREATION 

MECHANISM 


We have made attempts in this section to interpret the mechanism of particle creation as the 
phenomenon of Hawking like radiation from the homogeneous and isotropic FLRW space-time model. 
Normally, in the particle creation mechanism, the dissipative term behaves as an effective bulk viscous 
pressure, and, as a result, there is a negative pressure term in the Einstein’s field equations. Also, it 
is found that, by proper choice of the particle creation rate, there is an accelerated expansion of the 
universe both at the early stage (inflation), and at late-time, i.e., a complete description of the cosmic 
history. So, it is very natural to enquire whether there is any similarity between this mechanism with 
Hawking radiation as the inflationary stage can be described by the Hawking radiation in the FLRW 
universe 


39|, 


40( 


In Hawking radiation, initially due to the enormous size of the black hole, the evaporation process 
was very slow, and the process gradually became faster and faster with the diminision of the size of the 
black hole. As a result, the temperature of the black hole also increases with the process. At the end, 
when the black hole becomes of Planck dimension, quantum gravity should come into picture. On the 
other hand, the evolution of the universe is in the reverse direction of the black hole evaporation. At 
the beginning, the quantum gravity effects are important due to Planck size of the universe. But, with 
the evolution of the universe. Hawking radiation comes into picture, and the temperature gradually 
decreases. 

Moreover, there is another basic difference between these two evolution processes. In black hole 
evaporation, the created particles escape outside the event horizon, and, move to asymptotic infinity, 
while for evolution of the FLRW, the situation is just reversed — the particles created near the 
(apparent) horizon will move inside. As a result, due to black hole evaporation, there is a loss of 
energy, but the universe gains energy due to particle creations. Furthermore, due to isotropic nature 
of the FLRW model, the radiation should be uniform in all direction, and, we have from the Stefan- 


Boltzmann radiation law (SBRL) 


p _ 

dt 


<tAhT\ 


(39) 
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where a = j60h ^is the Stefan-Boltzmann constant, T is the radiation temperature, Q is 

the heat radiated by the black body, P is the net radiated power, and Ah is the radiating area. 


Using the above SBRL in the hrst law of thermodynamics, i.e.. 


we obtain 


dQ d dV 

^ = IT = 


(40) 


p + 3H{{p + p)-aT^} = 0. 


( 41 ) 


So, comparing with matter 
isentropic (i.e., using Eq. (jSjl), 


conservation equation ([5]), and considering the thermal process to be 
the particle creation rate is related to the temperature as 


r = 


7 ^’ 


Now, choosing E for the present model (in Eq. (117^ 1. we get 


(42) 


aT^ = + 3H^+ a‘^). ( 43 ) 

In the early phases of the evolution of the universe, H is very large, so, the above equation can 
approximately be written as 


aT ^ ~ 87/7^ = Kjp , 


(44) 


which is the usual black body radiation. On the other hand, at late-times, when the universe is 
very nearly to the equilibrium configuration, the Clausius relation becomes 


TS = Q = constant. (45) 

But, at late-time, the entropy is proportional to ‘a’ (see Eqns. (I20l) and ([22])), and, hence, S (x Hq a. 
Thus, from the above relation, we get 
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T oc 1/a, (46) 

which is nothing but the present CMB temperature. Lastly, it should be noted that, in the early 
phase of the universe, when F ~ Lf, T oc so, the temperature is not exactly Hawking temperature 

(r oc H), rather, Hawking-type radiation. The same is also true for the late-time evolution (where 
CMB temperature is dominant over Hawking temperature). Finally, the model becomes very close to 
the ACDM era; H is found to be a constant, and, hence, the particle creation rate, as well as, the 
temperature are constant, and there is no analogy with Hawking type radiation. 


VII. SUMMARY AND DISCUSSIONS 


Relativistic cosmology with usual perfect fluid (having barotropic equation of state) as cosmic 
substratum is considered in the present work in the context of non-equilibrium thermodynamics. In 
the framework of particle creation mechanism, dissipative phenomenon is reflected as an effective bulk 
viscous pressure. In the second order formulation of non-equilibrium thermodynamics due to Israel 
and Stewart, the dissipative pressure acts as a dynamical variable whose evolution is characterized 
by non-linear inhomogeneous evolution, however, the entropy flow vector satisfies the second law 
of thermodynamics. Due to our inability in solving the non-linear evolution equation, the thermal 
process is assumed to be adiabatic (i.e., the entropy per particle is constant), and, as a result, the 


dissipative pressure is linearly related with the particle creation rate. In earlier works 


23-1291], the 


particle creation rate at different cosmic stages had been chosen phenomenologically (with some basis 
from thermodynamics), and the solutions for different physical and thermodynamical variables show 


a continuity across the transition epochs. Subsequently, a single particle creation rate 


3ll | describes 


the whole evolution, but no analytic solution was found. The present work also considers the particle 
creation rate as a function of the Hubble parameter in a phenomenological way which not only describes 
the cosmic story from inflation to late-time acceleration, but also, the model has an analytic solution. 
Also, graphically, we have shown the evolution of the scale factor, Hubble parameter, deceleration 
parameter, the thermodynamic variables, namely, the temperature, entropy, and the number density. 
It is found that, at late-time, the model asymptotically approaches to ACDM, and it agrees with 
the latest Planck data set 3^. A cosmographic analysis has also been done for our model, and, the 


cosmographic parameters, namely, the jerk (j), snap (s), lerk (1), and the m parameter have been 
graphically shown against the cosmic time. From the figures, we conclude that the above parameters 
respectively have one, three, two, and two transitions through the cosmic evolution. Further, from the 







19 


field theoretic view point, it has been shown that, it is possible to have a minimally coupled scalar field 
as equivalent to the cosmic fluid for the description of the evolution of the universe. In the late-time 
asymptotic limit, the self similar potential behaves as cosmological constant. Finally, in resemblance 
with Hawking radiation, it is found that, at early epoch of the evolution, the temperature corresponds 
to black body radiation, while for very close to ACDM, the temperature is related to the CMBR — 
in both cases, the temperature is not the Hawking temperature, rather, Hawking type radiation. 
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